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Abstract 

We report an accurate study of interactions between Benzene molecules using variational quan¬ 
tum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these 
results with density functional theory (DFT) using different van der Waals (vdW) functionals. 
In our QMC calculations, we use accurate correlated trial wave functions including three-body 
Jastrow factors, and backflow transformations. We consider two benzene molecules in the parallel 
displaced (PD) geometry, and find that by highly optimizing the wave function and introducing 
more dynamical correlation into the wave function, we compute the weak chemical binding energy 
between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and 
-2.7(3) kcal/mol, respectively. The best estimate of the CCSD(T)/CBS limit is -2.65(2) kcal/mol 
[E. Miliordos et al, J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods 
give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to 
results from the best quantum chemistry methods. 
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I. INTRODUCTION 


Weak noncovalent van der Waals (vdW) interactions are fundamental to a wide range 
of topics relevance to physics, chemistry, and biology. A prototype vdW interaction is that 
resulting from the stacking between aromatic ringsi. This interaction is crucially important 
in biological systems such as protein folding^, DNA’s structure and stability^. In addi¬ 
tion, aromatic rings interactions play key roles in drug design^, electronics^, optical prop¬ 
erties of materials^, polymer stability^, conjugated carbon networks^, and crystal growth 
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processe9^“. 

In general vdW interactions are difficult to model accuratelyl^il^. Local and semilocal 
density functionals are unable to describe the long-range electronic correlation energy which 
is the main part of the vdW forces. We use quantum Monte Carlo (QMC) and modern 
non-local exchange-correlation (XC) functionals. Previously empirical and semi-empirical 
corrections were used. For instance, one approach was to add empirical, pairwise atomic 
dispersion corrections of the form -Cg/R®. To avoid double-counting electron correlation 
effects at short range, these contributions can be damped for small inter-atomic distances 
R. This method is referred to as DFT plus dispersion (DFT-D), and has been applied on 
different systems using various exchange-correlation (XC) funtionalsl^^— . 

Less empirical approaches were also developed. Effective nonlocal potentials were 
introduced^, where the parameters were determined using htting to ab initio results. Becke 
and Johnson obtained dispersion coefficients Cq,Cs, Cio from the multipole moments^S. 
In their work, the moments were obtained from an electron and its exchange hole. The 
reliability of these approaches is similar to the DFT-D formalism. 

Another density functional based method, widely applied on noncovalent systems, is the 
combination of DFT with symmetry-adapted perturbation theory, refereed as DFT-SAPT or 
SAPT(DFT)^^— . In these approaches, the dispersion term is obtained using the frequency 
dependent density susceptibility function of time dependent DFT (TD-DFT). The perturba¬ 
tion theory of intermolecular interactions can accurately predict the complete intermolecular 
potential energy surfaces for weakly bound molecular complexes^. From the point of view 
of perturbation theory, all the intermolecular interactions, including van der Waals interac¬ 
tions, contain four fundamental physical contributions: electrostatic, induction, dispersion, 
and exchange. The interactions differ only by proportions of these ingredients. The strongest 
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of those interactions involve a larger negative contribution from the electrostatic forces as 
in hydrogen bonded systems. If the electrostatic contribution is small, like in interactions of 
rare gas atoms with molecules, the minima depths are often below 1 kcal/mol. Therefore, 
from this point of view, a system such as the benzene dimer is also a van der Waals complex. 

Finally, by including nonlocal terms in DFT correlation energy functional, vdW-DF92^>2^ 
include the long range nonlocal correlation energy obtained by the plasmon pole approxima¬ 
tion. These functionals were originally applied on different van der Waals systems to obtain 
potential energy curve (PEC)^. DFT-vdW functionals result in signihcant improvements in 
equilibrium spacings between noncovalently bound complexes, as well as in binding energy 
of weak interacting systems. The efficiency and accuracy of different DFT-vdW functionals 
on bulk systems were analyzed^^— . The functionals are a clear improvement over semi-local 
functionals, although tests on a wider range of systems are desirable. 

Numerous theoretical works have compared different approximate quantum chemistry 
based methods for noncovalent weakly bound systems^^— . Particularly the coupled-cluster 
theory through perturbative triplets CCSD(T) which is often considered as the gold standard 
for chemical accuracy^*^^. However, due to its substantial computational cost, scaling as 
where N is the number of electrons, more efficient methods for vdW systems are highly 
desirable. These methods also cannot be applied to condensed matter. Thus we consider 
Quantum Monte Carlo (QMC), and compare with non-local density functionals. 

Quantum Monte Carlo, which solves the electronic Schrodinger equation stochasticallj^^— , 
is an alternative approach to quantum mechanical methods. Diffusion quantum Monte Carlo 
(DMC) provides accurate energies for vdW systema^^— . DMC is also able to produce an 
accurate description of systems where many-body interactions play a key role^i^. In general 
QMC based methods are faster than the most accurate post-Hartree-Fock schemes for large 
number of particles N. The computational cost of QMC methods scales usually as 
depending on the method. 

The benzene dimer has become a benchmark system for electronic structure methods 
for systems where van der Waals interactions are important. Despite its simplicity the 
problem of identifying the global minimum structure is particularly challenging as there are 
only subtle differences in the binding energies of the different conhgurations. According 
to quantum chemistry results, two critical factors for the binding energy of the benzene 
dimer are basis set and electron correlation:^. Our previous comprehensive study of benzene 
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molecules^ illustrates the importance of basis set in QMC energy calculations. Once the 
Jastrow factor is optimized by keeping fixed the Slater determinant, we obtained a good 
description of the atomization energy of the benzene molecule only when the basis of atomic 
orbitals is large enough and close to the CBS limit. In this work, we demonstrate that by 
using better trial wave functions and converged basis sets, we obtain a chemically accurate 
description of binding energy between aromatic rings. 


II. COMPUTATIONAL DETAILS 

We used the CASINO QMC code^ with a trial function of Slater-Jastrow (SJ) form, 

d^r(R.) = exp[J(R)] det[?/^n(rj)] det[?/>n(rj)], (1) 

where R is a 3A-dimensional vector that defines the positions of all N electrons, r ■ is the po¬ 
sition of the i’th spin-up electron, rj is the position of the j’th spin-down electron, exp[J(R)] 
is the Jastrow factor, and det['0„(rj)] and det['0n(i’j)] are Slater determinants of spin-up and 
spin-down one-electron orbitals. These orbitals were obtained from DFT calculations using 
the plane-wave-based Quantum Espresso code^. We used the local density approximation 
(LDA) to generate the orbitals in the Slater determinant for the trial wave function. We 
chose a very large basis-set cut-off of 200 Ry to guarantee convergence to the complete basis 
set limit^. The plane-wave orbitals were transformed into a blip polynomial basis^*^. The 
quality of the blip expansion , meaning the fineness of the blip grid, can be improved by 
increasing the grid multiplicity parameter and consequently results in a greater number of 
blip coefficients. The value of this parameter in our work is 2.0. The local density approxi¬ 
mation (LDA) pseudopotentials are generated using the OPIUM pseudopotential generation 
prograiii^. We also checked that the Kleinman-Bylander— transformation did not generate 
ghost states. In our DMC calculations, the pseudopotential energy was evaluated using a 
variational technique^. We used DMC time steps of 0.01 a.u. and 0.04 a.u. and extrapolated 
the results linearly to zero time step. 

The Jastrow factor is a positive, symmetric, explicit function of interparticle distances. 
We used a Jastrow factor consisting of polynomial one-body electron-nucleus (IB), two-body 
electron-electron (2B), and isotropic three-body electron-electron-nucleus (3B) terms. The 
main approximation in fermionic QMC is the fixed node approximation. To reduce this 


4 


error, we used backflow transformation (BF) in our trial wave functions^. In the backflow 
transformation, the orbitals in the Slater determinant are evaluated not at the actual electron 
positions, but a quasi-electron positions that are functions of all the particle coordinates. 
The backflow function, which describes the offset of the quasi-electron coordinates relative 
to the actual coordinates, contains free parameters to be determined by an optimization 
method. It allows the nodal surfaces to move within variational optimization, so with BF 
the QMC is no longer strictly fixed node. However, the subsequent DMC computations use 
the nodal surface that were determined during the VMC step. 

In QMC calculations, correlated wave function can be obtained by replacing the single 
determinant by a sum over configuration state functions (CSFs), using BF transformations 
of the electronic coordinates or by using pairing wave functions^^— . Our BF transformation 
includes both electron-electron and electron-proton terms, given by 

V({ri}) = r, + {''-'’({rj}) + (2) 

where Xj({rj}) is the coordinate of electron i which depends on the configuration of the 
system {rj}, and are electron-electron and electron-proton backflow 

displacements of electron i, respectively, given by 

Afe 

Np 

(4) 

I 

where aijivij) and (^uiru) are polynomial functions of electron-electron and electron-proton 
distance, respectively, containing optimizable parameters. 

We use two methods for wave function optimization: variance minimization and energy 
minimizatiort^*^. The parameters of Jastrow and backflow are first optimized by variance 
minimization at the variational Monte Carlo (VMC) level^^— . Since trial wave functions 
generally cannot exactly represent an eigenstate, the energy and variance minima do not co¬ 
incide. Therefore energy minimization should produce lower VMC energies. We have found 
that the lower VMC energies lead to lower DMC energies if we use backflow transformations. 
This is due to improved many-body nodes as well as reduction in the errors induced by using 
non-local pseudopotentials. 
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FIG. 1. (colour online) Parallel-displacement (PD) geometry of benzene dimer {C 2 h symmetry) 
studied in this work. R indicates the distance between centers R = y^(Rl'l-^|i. In our study 
centers displacement Ri = 1.6 A is fixed. 

We used the Quantum Espresso code^ for DFT-vdW calculations with ultrasoft pseudopotentials^ 
and Perdew-Burke-Ernzerhof (PBE)^ exchange correlation functionals. The plane wave ba¬ 
sis had a well-converged cut-off of 80 Ry. 


III. RESULTS AND DISCUSSION 

We used the experimental geometry for the benzene molecule, where the C — C and 
C — H bond lengths are 1.39 and 1.09 A, respectively. Experiments support the existence 
of these three benzene dimer conhgurations^, the parallel (sandwich), the T-shaped (C 2 w), 
and the slipped-parallel or parallel-displacement (PD, C 2 h}^ conhgurations. The T-shaped, 
where two benzene molecules are perpendicular to each other, and the PD conhgurations 
are more energetically favored than the parallel sandwich geometry^^^— . T-shaped and PD 
conhgurations of benzene dimer are almost isoenergetic, and the benzene dimer potential 
energy surface is quite hat with several local minima separated by tiny barriers. In this work 
we focus on the PD conhguration as shown in Figure [T] Centers of two parallel benzene 
rings are displaced by Ri = 1.6 A and is hxed in our calculations. 
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Figure [2] illustrates DFT potential energy curve obtained using different vdW function¬ 
als. We use vdW-DFl^, vdW-DF3^, vdW-DF-obk8, vdW-DF-ob86, vdW-DF2-B86R-^‘^, 
vdW-DF-C09, vdW-DF2-C09^, and vdW-DF-cx— functionals. All vdW functionals use 
Slater exchange and PW— correlation functionals. The non-local terms are either vdW-DFl 
or vdW-DF2. Employing various gradient correction on exchange energy is the main dif¬ 
ference between these functionals. Using polynomial htting (Appendix), the optimal DFT 
binding energies obtained by vdW-DFl, vdW-DF2, vdW-DF-obk8, vdW-DF-ob86, vdW- 
DF2-B86R, vdW-DF-C09, vdW-DF2-C09, and vdW-DF-cx are -3.1, -2.8, -3.1, -3.2, -2.4, 
-3.0, -1.5, -2.9 kcal/mol at R = 3.7, 3.65, 3.58, 3.60, 3.63, 3.57, 3.71, and 3.65 A, respec¬ 
tively. Free-energy landscape calculations using Car-Parrinello molecular meta-dynamics 
methods using the BLYP density functional with dispersion corrections predict that T- 
shape geometry is more stable at all temperatures^. However, the PD conhguration with 
C 2 h symmetry has been determined using optical absorption spectroscopy, whereas a polar 
V-shape conhguration with C 2 ,; symmetry has been suggested by multiphoton ionization 
mass spectroscopy^!^. 

TABLE I. QMC energies (kcal/mol) of the PD benzene dimer obtained by VMC and DMC using 
one-body (IB), two-body (2B), three-body (3B) Jastrow factors and Backflow(BF) transformations. 
Energies are calculated at different distance geometry R A. Energy differences are calculated with 
respect to the large R = 10 A. Experimental binding energy is -2.4(4) kcal/mol^. 


R VMC DMC 

VMC 

DMC 

VMC 

DMC 

m+2B lB-h2B 

lB-h2B m+2B 

1B+2B 

1B+2B 


-h3B 

-F3B 

-h3B-hBF -h3B-hBF 

3.0 7.23(9) 5.63(9) 

6.9(2) 

5.3(1) 

5.0(2) 

4.4(3) 

3.5 2.50(9) -0.81(9) 

0.0(2) 

-0.9(1) 

-1.2(2) 

-1.8(2) 

4.0 0.55(9) -1.45(9) 

-0.6(2) 

-1.5(1) 

-2.2(2) 

-2.5(2) 

5.0 0.02(9) -0.71(8) 

-0.0(2) 

-0.8(1) 

-0.6(2) 

-1.1(2) 

6.0 0.00(9) -0.31(9) 

0.0(2) 

-0.4(2) 

-0.3(3) 

-0.4(2) 


Table [H lists QMC energies of benzene dimer in the PD conhguration obtained by VMC 
and DMC methods at diherent separation distance R. Lowest VMC and DMC energies are 
obtained at R = 4.0 A. Adding 3B-Jastrow factor substantially improve the VMC energies. 
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FIG. 2. (colour online) DFT energy of PD benzene dimer as a function of their separation obtained 
using different vdW-DF functionals. The reference is chosen at R = 10 A. Using polynomial fitting 
(Appendix), the optimal DFT binding energies obtained by vdW-DFl— , vdW-DF22S, vdW-DF- 
obk8, vdW-DF-ob86, vdW-DF2-B86R-^>^, vdW-DF-C09, vdW-DF2-C092^, and vdW-DF-cx^ are 
-3.1, -2.8, -3.1, -3.2, -2.4, -3.0, -1.5, -2.9 kcal/mol at R = 3.7, 3.65, 3.58, 3.60, 3.63, 3.57, 3.71, and 
3.65 A, respectively. The experimental binding energy is -2.4(4) kcal/mol^^. 


The 3B-Jastrow function takes care of what is missing in the IB and 2B Jastrow factors, 
meaning, the explicit dependence of the electron correlation on the ionic positions. At 
R = 4.0 A, the difference between VMC energies obtained with and without 3B-Jastrow 
factor is about -1.1(2) kcal/mol, wheras in DMC the difference is negligible. However, BF 
transformations significantly lower the QMC energies. At R = 4.0 A, BF transformations 
lower the energies by -1.6 and -1.0 kcal/mol at VMC and DMC calculations, respectively. 
Adding 3B-Jastrow factor and BF transformations improves the VMC and DMC energies 
by about -2.7(2) and -1.0(2) kcal/mol, respectively. 













































FIG. 3. (colour online) QMC energy of benzene dimer as a function of their distance ob¬ 
tained using VMC and DMC methods. Energy differences are calculated with respect to the 
large R = 10 A. Using polynomial fitting (Appendix), the optimal values of binding energies 
obtained by DMC/1B+2B, VMC/1B+2B+3B, DMC/1B-|-2B-|-3B, VMC/1B-|-2B-|-3B-|-BF, and 
DMC/1B+2B+3B-I-BF are -1.7(2), -0.9(2), -1.8(2), -2.3(4), and -2.7(3) , respectively. 

The VMC energies obtained with BF are lower than DMC energies without BF. This 
indicates that VMC with BF could be a useful level of theory to describe nonlocal long 
term interactions. In general, VMC calculations are significantly less expensive than DMC 
ones. Also, VMC has advantages for calculating expectation values of quantities more than 
the energy. Our QMC results clearly demonstrate that increasing the complexity of the 
wavefunction by including the BF correlations and 3B-Jastrow terms, plays key role in 
accurate describing of vdW interactions. 

Figure E] illustrates QMC potential energy curve of benzene dimer as a function of their 
distance obtained by VMC and DMC calculations. Using polynomial fitting (Appendix), 
the optimal values of binding energies at VMC level obtained by lB-|-2B-|-3B-Jastrow and 
lB-|-2B-|-3B-Jastrow plus BF correlations are -0.9(2), and -2.3(4) kcal/mol, respectively. 
DMC optimal values obtained by lB-|-2B-Jastrow, lB-|-2B-|-3B-Jastrow, and 1B-I-2B-I-3B- 
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Jastrow plus BF are -1.7(2), -1.8(2), and -2.7(3) kcal/mol, respectively. Using only 1B-I-2B- 
Jastrow factor, VMC is unable to provide a bound benzene dimer. Using the same Jastrow 
factor the DMC binding energy is close to those ones obtained by S0S-MP2 method^. 
Employing 3B-Jastrow factor signihcantly improves the VMC binding energy, whereas it 
doesn’t lower the DMC binding energy considerably. At the variational level, the inclusion 
of a 3B-Jastrow term provides additional dynamical correlation into the wavefunction and 
it is essentially useful for studying nonlocal vdW interactions. 

Although a 3B-Jastrow factor improves the binding energy of benzene dimer in VMC, 
substantial enhancement is obtained by employing BF correlations. The optimal values 
indicate that VMC and DMC energies are improved by -1.4 and -0.9 kcal/mol, respectively. 
It suggests that BF is effective at improving the nodal surface of benzene dimer in PD 
geometry. Considering that the LDA wave function often has too many nodal pockets, it is 
conceivable that BF coordinate transformations could modify the number of nodal pockets 
of a wave function. 

Table m shows the binding energy of benzene dimer in PD geometry obtained by dif¬ 
ferent methods. Since the benzene dimer is a standard test for high-level quantum chem¬ 
istry methods for proper characterization of vdW interactions there are many more results; 
a comprehensive comparison between different high-level ab initio approaches is reported 
recently^. They provide databases for noncovalent interactions. 49 bimolecular complexes 
in 345 geometry conhgurations are partitioned into subsets based on bonding motif. Benzene 
dimer is in the dispersion-dominated subset. Our DMC result is close to CCSD(T) (Table 
ED, which is considered the gold standard for chemical accuracy. The complete basis set 
(CBS) limit can now be estimated more precisely in the CCSD(T) framework. Our VMC re¬ 
sults is comparable with those ones obtained by LRDMC method using Jastrow-AGP wave 
function^. The correlated antisymmetrized geminal power (AGP)^ is the particle num¬ 
ber conserving version of the Bardeen-Cooper-Schrieffer (BCS) wave function. A singlet 
valence bond between two electrons of opposite spin is determined by a geminal function. 
This framework successfully applied to identify the Kekule and Dewar contributions to the 
chemical bond of the benzene molecule^. 

The T-shape conhguration was studied before using QMC methods^. Using IB, 2B, and 
3B Jastrow factors, they found that the binding energy of T-shape conhguration obtained 
by FN-DMC is -2.88(16) kcal/mol. An accurate study of noncovalent systems illustrates the 
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TABLE II. Binding energies of the PD benzene dimer obtained at different level of theories. 
The zero point vibrational energy is not included. 


Method Basis set 

R (A) i7;,(kcal/mc 

»1) References 

VMC/lB-t2B-t3B-tBF 

3.9(3) 

-2.3(4) 

This work 

DMC/lB-t2B-t3B-tBF 

3.8(3) 

-2.7(3) 

This work 

DFT-D/BLYP|TZVP 

3.486 

-2.88 

Ref.^ 

vdW-DF2 

3.65 

-2.8 

This work 

vdW-DFl 

3.7 

-3.1 

This work 

vdW-DF-obk8 

3.58 

-3.1 

This work 

vdW-DF-ob86 

3.60 

-3.2 

This work 

vdW-DF2-B86R 

3.63 

-2.4 

This work 

vdW-DF-C09 

3.57 

-3.0 

This work 

vdW-DF2-C09 

3.71 

-1.5 

This work 

vdW-DFT-cx 

3.65 

-2.9 

This work 

CCSD(T)|CBS 

3.9 

-2.65(2) 

Ref. 49 

CCSD(T)|CBS(AaDZ) 

N/A 

-2.73 

Ref. 32 

MP2|CBS 

3.66 

-5.00(1) 

Ref. 49 

JAGP-LRDMC 

4.1(2) 

-2.2(3) 

Ref. 44 

FNDMC 

N/A 

-1.65(42) 

Ref. 45 

Experiment 

N/A 

-2.4(4) 

Ref.^ 

Experiment 

N/A 

-1.6(2) 

Ref.^ 


importance of Jastrow factor optimization in obtaining reliable results^. They investigated 
in detail all technical parameters of QMC simulations. They also have found that the 
binding energies for T-shape and PD configurations are -3.77(39) and -1.65(42) kcal/mol, 
respectively. The CCSD(T) estimates of T-shape binding energy is -2.74 kcal/mol^. By 
considering the reduction of the binding energy due to the ZPE (AZPE=0.37 kcal/mol)^, 
our DMC energy obtained by 3B-Jastrow and BE correlations is in excellent agreement with 
experiment^. Among DFT results obtained by different vdW functionals, vdW-DF2-B86R 
energy is close to our VMC/1B-|-2B-|-3B-|-BF result. vdW-DFT-cx and vdW-DF2 energies 
are close to our DMC/1B-I-2B-I-3B-I-BF energy. 
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We found that the dependence of the DMC energies on the quality of the trial wave 
function is significant. Whereas the Jastrow factors keep electrons away from each other 
and essentially improve the trial wave function, they do not change the nodal surfaces. It 
has been argued that BF transformation and 3B Jastrow correlation arises as the next-order 
improvements to the standard Slater-Jastrow wave function^. As our results show, the 
importance of BF correlations within DMC calculations is that they alter the nodal surface 
and can therefore be used to reduce the FN error. However, more complexity of BF-WF 
comparing to SJ-WF causes additional computational cost in QMC calculations. One of 
the most expensive operations in QMC calculations is the evaluation of the orbitals and 
their first two derivatives. The evaluation of the collective coordinates in BF-WF introduce 
significantly more computational cost. Moreover, whereas QMC calculations with SJ-WF 
require only the value, gradient, and Laplacian of each orbital '0, BF calculations also require 
cross derivatives such as /dxdy. The most important complicating factor arising from 
BF transformations is that they make each orbital in the Slater determinants depend on the 
coordinates of every particle. 

Figure |4] illustrates the potential energy curves for PD benzene dimer as function of 
ring centers distance at different level of theory. CCSD(T) and MP2 results were calculated 
using aug-cc-pVQZ* (=aug-cc-pVQZ - g function on carbon - f function on hydrogen) basis 
set^. At the MP2 level, the equilibrium distances between aromatic rings are Ri = 1.6 and 
R 2 = 3.4 A. At the CCSD(T) level, they are Ri = 1.6 and R 2 = 3.6 A. MP2 calculations 
produce over-binding. Comparing to T-shaped benzene dimer conhguration, which we have 
not studied in this work, it has been claimed that PD configuration has larger electrostatic 
interactions^^, since the positive hydrogens on each aromatic ring are located on top of the 
negative carbons of the other ring. In the PD geometry two benzene rings are closer together 
than in either the sandwich or the T-shaped configurations. Therefore, the incursion of the 
electronic orbitals of each ring make the electrostatic interaction more stable. 


IV. CONCLUSION 

We hnd that QMC can give chemical accuracy for the benzene dimer. The good agreement 
among our results, experiments and quantum chemistry methods, is an important sign of 
the capability of the QMC based methods to provide accurate description of very weak 
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FIG. 4. (colour online) Energy for PD benzene dimer as function of their centers dis¬ 
tance obtained by different methods. We only compare our VMC/1B-|-2B-|-3B-|-BF (VMC), 
DMC/1B-I-2B-I-3B-I-BF (DMC) results. The MP2, and CCSD(T) data were taken from Ref 92. 


intermolecular interactions based on vdW dispersive forces. We find that adding 3B-Jastrow 
terms and BF transformations leads to significant improvement in the accuracy of the weak 
vdW interaction between aromatic rings. BF-VMC energies are significantly lower than 
SJ-VMC and therefore BF-VMC could be useful alternative for a SJ-DMC calculations 
, which are more expensive. The accuracy of our VMC results compared with DMC is 
evidence of the high accuracy of our trial wave function. BF correlations give substantial 
enhancement in trial wave function of aromatic rings. By improving the nodal surface of 
wave function, it leads to a signihcant reduction in binding energy between two benzene 
molecules. Improved trial wave functions will be useful in VMC calculations of quantities 
other than the energy, which are usually more difficult to obtain accurately than the energy. 
We used single determinant wave functions in this work , but BF can also be combined with 
other types of wave functions such as multideterminant or pairing wave functions. 
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V. APPENDIX: ENERGY CURVE FITTING FUNCTIONS 


In this appendix we report fitting parameters. We used piecewise polynomial fitting 
functions. Tables uni and nYi list fitting parameters including smoothing parameter a, 
the sum of squares due to error (SSE), R-square, Root Mean Squared Error (RMSE), and 
optimal parameters {opt) for each method. 

TABLE III. Piecewise polynomial fitting parameters for QMC results. Smoothing parameter a, 
the sum of squares due to error (SSE), R-square, Root Mean Squared Error (RMSE), and optimal 
parameters (opt) for each method. Energies E and distance geometry R are in kcal/mol and A, 
respectively. 


VMC 

DMC 

VMC 

DMC 

VMC 

DMC 

1B+2B 

1B+2B 

m+2B 

1B+2B 

1B+2B 

1B+2B 



+3B 

+3B 

-h3B-hBF 

+3B+BF 


a 0.9989716 0.9999155 0.9999155 0.9999155 0.9998607 0.9998607 

SSE 0.00372 0.0002223 0.0002638 0.0002345 0.0003988 0.0005139 


R-square 

0.9999 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

RMSE 

0.1694 

0.1379 

0.1502 

0.1416 

0.1443 

0.1638 

^opt 

0.00 

-1.66 

-0.87 

-1.77 

-2.31 

-2.71 

^opt 

4.50 

3.78 

3.78 

3.80 

3.85 

3.79 
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TABLE IV. Piecewise polynomial fitting parameters for DFT results. Smoothing parameter a, 
the sum of squares due to error (SSE), R-square, Root Mean Squared Error (RMSE), and optimal 
parameters (opt) for each method. Energies E and distance geometry R are in kcal/mol and A, 
respectively. 



vdW- 

vdW- 

vdW- 

vdW- 

vdW- 

vdW- 

vdW- 

vdW- 


DF- 

DF- 

DF2- 

DF- 

DF2- 

DF- 

DEI 

DF2 


obk8 

ob86 

B86R 

COO 

COO 

cx 



a 

0.9996214 0.9996214 0.9996214 0.9996214 0.9996214 0.9997703 0.99977035 0.99977035 

SSE 

0.001203 

0.001034 

0.0011 

0.0007068 0.0008411 0.0003853 

0.0007318 

0.0007439 

R-square 

0.9998 

0.9998 

0.9997 

0.9998 

0.9998 

0.9999 

0.9999 

0.9999 

RMSE 

0.1539 

0.1426 

0.1471 

0.1179 

0.1286 

0.111 

0.1529 

0.1542 

^opt 

-3.13 

-3.16 

-2.37 

-3.03 

-1.46 

-2.90 

-3.06 

-2.79 


3.58 

3.60 

3.63 

3.57 

3.71 

3.65 

3.71 

3.65 
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